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Abstract 

We use previously published high-resolution synchrotron polarization data to perform 
an angular dispersion analysis with the aim of charactering magnetized turbulence in 
M51. We first analyze three distinct regions (the center of the galaxy, and the northwest 
and southwest spiral arms) and can clearly discern the turbulent correlation length scale 
from the width of the magnetized turbulent correlation function for two regions and 
detect the imprint of anisotropy in the turbulence for all three. Furthermore, analyzing 
the galaxy as a whole allows us to determine a two-dimensional Gaussian model for the 
magnetized turbulence in M51. We measure the turbulent correlation scales parallel and 
perpendicular to the local mean magnetic field to be, respectively, 5|| = 98 it 5 pc and 
5_L = 54 lb 3 pc, while the turbulent to ordered magnetic field strength ratio is found 



to be Bt/Bo = 1.01 it 0.04. These results are consistent with those of iFletcher et al. 



(120 111 ), who performed a Faraday rotation dispersion analysis of the same data, and our 



detection of anisotropy is consistent with current magnetized turbulence theories. 



Subject headings: galaxies: spiral - galaxies: magnetic fields - galaxies: ISM - galaxies: 
individual: M51 
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Introduction 



The magnetized diffuse interstellar medium, of the Milky Way and other galaxies, is turbulent 
and so an understanding of its properties and role should include the quantities commonly used to 
describe turbulence: characteristic length scales, power spectra, the relative energies in the mean 
and fluctuating components, and so on. One important property of magnetohydrodynamic (MHD) 
turbulence is that the random fluctuations in the inertial range are not necessarily isotropic, as is 
the case in the classical picture of purely hydrodynamic, incompressible Kolmogorov turbulence: the 
correlation length of magnetic fluctuations can be larger along the mean field direction compared 
to the perpendicular direction. This mean field can be either an external large-scale magnetic field 
or simply the magnetic field at the largest scale of a turbulent eddy acting on fluctuations within 
the eddy on smaller scales. As well as the inherent anisotropy of MHD turbulence dynamical effects 
in the ISM flow such as shocks and shear, due to localized sources such as supernovae or global 
features like differential rotation, can also imprint anisotropy on the turbulence. 

There have been a few indications f rom o bservations that magnetic field fluctuations in the ISM 
exhibit anisotropics. iBrown fc Taylod (l200lh binned Faraday rotation measures (RM) for extra- 
galactic (EG) sources in the Galactic plane in the range 82deg < I < 146 deg and found that the 
variance in a bin is correlated with the magnitude of the mean RM; higher RMs are associated with 
stronger fluctuations, and it was proposed th at this occurs bec ause the fluctuations in the magnetic 
field are mainly aligned with the mean field. iJaffe et al.l ()20ld ) fitted a model magnetic field to the 
synchrotron emission and EG RMs along the Galactic plane and found that an anisotropic (or in 
their terminology an ordered random) magnetic field component was required to fi t the observations 
along with both a mean field and an isotopic random magnetic field. Similarly, IJansson &: Farrai 
(|2012l ) required either an anisotropic (or striated field in their terminology) magnetic field com- 
ponent, or a correlation between the mean magnetic field and cosmic ray density, to obtain good 
fits to all-sky observations of synchrotron emission and RMs (in their model anisotropic random 
fields and close cosmic ray- t o-mea n magnetic field coupling are degenerate parameters). Away from 
the Milky Way, Beck et al.l (j2005l ) compared the observed increase in both total and polarized syn- 
chrotron emission at the strong shock fronts along the bars of the galaxies NGC 1097 and NGC 
1365 with theoretical expectations based on the compression and shear of random and mean mag- 
netic fields; the ir results indicate tha t strong anisotropic random magnetic fields are produced at 
these positions. iFletcher et al.l (j201ll ) attributed the order of magnitude difference between ordered 
magnetic field strengths obtained via equipartition estimates and Faraday rotation modeling, to a 
strong anisotropic random magnetic field in the nearby galaxy M51: this component is responsible 
for the strong polarized signal but contributes little to the Faraday rotation. 



In this paper we perform an angu lar dispersion analysis, ba sed on the work of lFalcet a- Gonial ves et al 



(|2008h . lHildebrand et all ^200% . and 



Houde et al 



(2009, 



2011 



on the Effelsberg 100-m/VLA A6.2 
cm synchrotron polarization map of lFletcher et al.l (j201ll ) (4" FWHM resolution and 1" sampling). 
We show in Figure [T] the global view of M51 in polarized flux provided by these data. There are 
clearly only three regions that can be used, or combined, for a dispersion analysis: the spiral arms 
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Fig. 1. — Polarized flux at A6.2 cm for M51. There are three regions that can be independently 
used (or combined) for a dispersion analysis: the spiral arms in the northeast and southwest, and 
the center of the galaxy. These regions are contained within the corresponding three red rectangles 
in the figure. The map is centered at RA (J2000) = 11^29"'52!4, Dec (J2000) = 47°11'43'.'5, and the 
contours are drawn at 20 to 80 percent (10 percent increments) of the peak polarized fiux density 
(173 /iJy/beam). 
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in the northeast and southwest, and the center of the galaxy. These regions are contained within 
the corresponding three red rectangles in the figure. 

Although these data were obtained with high spatial resoluti on, they will not al low us to get 



a handle on the magnetized turbulence power spectrum as done inlHoude et al.l (l201ll ) for Galactic 
molecular clouds (note that l" ~ 37 pc). Nonetheless, iFletcher et al.l (120 111 ) calculate from a 



Faraday rotation dispersion analysis that the size of a turbulence cell should be approximately 50 
pc. This suggests that we may be able to measure and determine this value independently with 
the data at hand, given the expected size of a cell and the aforementioned spatial resolution. We 
would be in a position to not only determine the number of turbulence cells contained in the average 
column of gas subtended by the telescope beam, but also get accurate values for the ratio of the 
turbulent-to-ordered magnetic energy in different parts of M51. Furthermore, in view of the large 
numbers of polarization measurements available with this map the statistics may be good enough 
to allow a study of possible anisotropy in the autocorrelation function of magnetized turbulence. 

We start in the next section by first giving a summary of the dispersion analysis in Section [21 a 
description of the data used for our analysis is given in Section [3l which is then followed in Section 
4.11 bv an isotropic dispe rsion analysis on the three a vailable regions in the manner presented in 



Houde et al. 



2011 



and iHildebrand et al.l (120091 ). A first attempt at measuring any potential 
anisotropy is presented in Section [4.2.11 through the independent analyses on displacement vectors 
that are grouped into two sets, which are either oriented approximately parallel or perpendicular 
to the local mean magnetic field. The derived turbulence autocorrelation functions can then be 
compared and any differences in their widths will reveal an anisotropy. Finally, in Section 14.2.2! 
we apply the dispersion analysis to M51 as a whole (i.e., by simultaneously using all available 
polarization vectors) to map the two-dimensional turbulence autocorrelation function to once again 
reveal any anisotropy in the magnetized turbulence, but in more detail. We finish with a discussion 
of our results in Section [5l a brief summary in Section [6l while more details concerning the dispersion 
analysis will be found in the Appendix. 



2. Angular Dispersion Analysis 



Structure functions have long been used in physics and astrophysics to characterize turbulence, 
as they allow for the treatment of power-law energy spectra, such as those found in Kolmogorov 



turbulence, without the mathematical divergen ces associated with stationary signal (IFrisch 



1995 



Beck et al .11 19991 : iFalcet a- Gonial ves et al.ll2008l ). Such structure functions, of varying orders, can 
be calculated for a range of physical parameters (e.g., velocity and density fields). In this pa- 
per, we intend to apply the angular dispersion analysis previously introduce d in the literature 
(jFalceta-Gongalves et aPbood ; lllildebrand et ahlbood : iHoude et aUbood boill ). where the chosen 
parameter is the orientation of the projection of the magnetic field on the plane of the sky. More 
precisely, we will use the polarization angle orientation in lieu of that of the magnetic field. For 
the polarization of synchrotron (or dust) emission this angle is orientated at 90° from that of the 
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projected magnetic field. Although we will provide a summary of the important relations required 
for the angular disp ersion analysis lat e r in S ection 12.21 a simplified exposition based on material 
that can be found in iHildebrand et al.l (|2009l ) is first given here. 



2.1. Angular Structure Function 

We start by defining the difference A$ {£) in the orientation of the magnetic field (unless 
otherwise specified, in this paper we will only concern ourselves with the plane of the sky component 
of the magnetic field) at two points separated by a distance i 



A$(€)EE$(r)-$(r + €), (1) 

with $ (r) the magnetic field orientation at position r (both r and £ are also understood to be 
located in the plane of the sky). Given a set of measurements <I> (r) on a polarization map, we can 
also define the following (second order) angular structure function 



(A$2 (£)) 



N{e) 



N{i) 
i=l 



(2) 



where (•••) denotes an average, N (i) is the number of pairs of field orientation measurements 
separated hy i = \i\, and stationarity and isotropy were assumed (i.e. , the structure function is only 
dependent on the magn itude of ^, and not on its orientation or r; see lFalceta-Gongalves et al.ll2008l : 
Hildebrand et ahlbood ;). 



The main assumption in our analysis consists in modeling the magnetic field, and therefore 
its orientation through $ (r), as being composed of a large-scale, ordered component $o (r) and a 
smaller scale, zero-mean, turbulent component $t (r). That is, we write 



$ (r) = (r) + cDo (r) , (3) 
which, if we further assume these two components to be statistically independent, leads to 



(A$2(^)> = (AcI>?(£)> + (Aci>2(£)). (4) 

We thus find that the structure function is composed of two angular components stemming from 
the contributions of the turbulent and ordered magnetic fields. Our assumption on the difference 
between the two scales therefore allows for their separation and analyses. This is exemplified in 
Figure [2] where we show (Panel a)) hypothetical turbulent and ordered contributions to the total 
angular structure function (Panel b), solid curve), as expressed in Equation (|4]). Panel c) presents 




Fig. 2. — An example of an angular structure function. Panel a) shows hypothetical turbulent 
(A$^ (^)) (triple-dot-broken curve) and ordered (A$q (dot-broken curve) contributions to the 
total angular structure function (A$^ (£)) plotted in Panel b) (solid curve). Panel c) presents the 
same information as Panel b) but displayed as a function of i"^ instead of i. The total angular 
structure function is subtracted from a Taylor series fit obtained from its values at i > 6" (dot- 
broken curve in Panels b) and c)) to get the turbulent autocorrelation function shown in Panel d). 
The correlation length characterizing turbulence is represented by S. 
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the same information as Panel b) for <^A$^ but displayed as a function of £^ instead of I. This 
is to show that the separation of the two length scales is sometimes easier to visualize, through the 
abrupt change in the slope of ^A<&-^ (^))) by using the square of the distance for the abscissa; we 
will use both representations for the data analyzed later in this paper. The total angular structure 
function (of Panels b) or c)) is the input to our problem as obtained from a polarization map, which 
we seek to analyze in order to characterize magnetized turbulence. 

The behavior of turbulent and ordered contributions to the total angular structure function of 
Panel a) can be qualitatively understood from the fact that, evidently, they must equal zero when 
£ = and then initially increase with i. The turbulent structure function will keep increasing until 
it reaches values for I that sufficiently exceed the turbulence correlation length 5, at which point it 
will reach its maximum value. This can also be understood more quantitatively with the following 
relation 

(A$?(£)) = ([$ar)-$t(r + ^)]') 

= 2[($2(0)>-(ci>2(^))], (5) 

where the (stationary and isotropic) turbulent autocorrelation function is defined with 

{^l{l))^{^,{r)^,{v + i)). (6) 

It therefore follows that (A<I>^ (£)) = 2 {^1 (0)) when £ ^> 5, as the turbulence is not correlated on 
such scales and its autocorrelation vanishes. The ordered structure function is expected to rise, at 
first monotonically, with increasing values of i in view of its larger-scale nature. For the example of 
Figure [2] we have characterized the turbulent component with a Gaussian autocorrelation function 
of width, or correlation length (i.e., its standard deviation) 5 = 2", while the ordered structure 
function was modeled with a Taylor expansion in powers of l"^ . This restriction to even powers in i 
is dictated from the assumption of isotropy for the structure function. 

With the previous assumption in the difference of the two length scales it becomes possible to 
model the ordered component (A<1>q (^)) independently of (A$^ (£)) by using values of i sufficiently 
large (i.e., sufficiently greater than 5) where any variation in (A<1>^ (£)) is negligible. For our 
example, we chose I > 6" to obtain the Taylor series fit given by the dot-broken curve in Panels b) 
and c) of Figure [2l This curve is then representative of (A<1>q (^)) but shifted up by the constant 
level of the turbulent component present in that range. More precisely if we define a function x (^) 
for this curve, then we write 



x(^) = 2(cI>2(0)> + (Aci>2(^)>. (7) 
In this paper we will focus on characterizing magnetized turbulence and we are thus interested in 
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isolating the turbulent component of the structure function, or, alternatively, its autocorrelation. 
The latter (multiplied by a factor of two) is readily evaluated from Equations dH) , ([5]) , and ([7]) with 



2($2(£)) = ^(£)_(AcI>2(£)> (8) 

and shown in Panel d) of Figure [2] (i.e., as the subtraction of the solid curve from the dot-broken 
curve in Panel b)). 

Although we could very well use the angular structure function in the manner presented in 
this section for our analyses of the M51 data, we will nonetheless for the rest of our discussion 
focus instead on the angular dispersion function 1 — (cos [A<I> (^)]). We note, however, that the 
properties, method, and technique discussed thus far for the structure function apply just as well to 
the dispersion function. In fact, we note that the two angular functions are simply related through 



l-{cos[A^{i)])c^^{A<^^{i)) (9) 

when A$ {£) <C 1. The advantage of the dispersion function is its close connection to the autocor- 
relation of the magnetic field (see Equation (llOp below) , which then naturally l eads to the study o f 



the magnetized turbulent power spectrum through a simple Fourier transform (iHoude et al.ll201ll ). 



Finally, it is important to note that in general the width of the turbulent autocorrelation 
function (and of the associated structure/dispersion function) is not solely due to the intrinsic 
correlation length of turbulence. A correlation scale brought about by the finite spatial resolution 
with which observations are realized will also combine to the intrinsic turbulent correlation length 
to set the overall width of the turbulent autocorrelation function. A further complication results 
from the related problem of signal integration through the line of sight. As we will soon see, 
a careful analysis of such effects will not only allow us to disentangle the intrinsic correlation 
length characterizing magnetized turbulence from the overall width of the turbulent autocorrelation 
function, but also to determine the level of turbulent energy contained in the medium probed by 
the observations. For this to be feasible, however, some approximation must be made on the nature 
of the turbulent autocorrelation function. T he case of iso t ropic Gaussian turbulence and beam 



profile functions was treated in the details in iHoude et al.l (|2009l ) (see their Section 2) and their 



main relations detailing the combination of the two length scales for the analysis of turbulence are 
given in the next section (see Equations (l22|) -( l24|l below). In this paper, we further provide an 
analysis of the more general case of anisotropic Gaussian turbulence in the Appendix, while the 
corresponding results are also summarized in Section 12.21 (see Equations (jl9p -( f2T]) ). This will in 
turn make possible the measurement of anisotropy in the turbulent autocorrelation function, which 
is another important parameter for the characterization of magnetized turbulence. 
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2.2. Angular Dispersion Function 

As previously mentioned, the analysis of the angular dispersion fu nction found i n the Appendix 



and summarized in this section follows that presented in Section 2 of iHoude et al.l (|2009l ) with the 
difference that we now allow for the presence of anisotropy in the turbulence. As stated in Section 
12.11 we are interested in the function cos [A<I> {£)] that is related to the magnetic field autocorrelation 
function through 

(cos[A^(£)]) = ; ) (10) 

^ ^ ^ (B-B(0)> 

where (• • • ) denotes an average and (B-B(£)) = (B (r) •B(r + £)). It is important to note that 
the magnetic field B is a weighted average (with the polarized fiux) through the thickness of the 
column of gas probed (i.e., the disk of M51) and across the telescope beam (see Equation (lAip ). 
The local, non-averaged, magnetic field B (x) at a point x is composed of an ordered field Bo(x) 
and a turbulent component Bt (x) such that 



B(x) = Bo(x)+Bt(x). (11) 



As was the case earlier, the displacement vector £ in Equation (llOp . and others that will 
follow, is understood to be located in the plane of the sky. We will further break £ down into two 
perpendicular components 



£ = £i+£2, (12) 

where, unless otherwise noted, £i and £2 are taken to be respectively perpendicular and parallel 
to the projection of the ordered component of the magnetic field Bq on the plane of the sky. For 
everything that follows statistical independence between Bo(x) and Bt (x), homogeneity in the 
strength of the magnetic fields (-Bq) = (Bo'Bq (0)) and (-B^) = (BfBt (0)), as well as, more 
generally, stationarity in the autocorrelation functions (Bo-Bo(£)) and (BfBt(£)) are assumed. 
The assumption of homogeneity in the field strength in particular, while clearly an idealization, is 
needed for securing a quantitative measure of turbulence from our data (see the Appendix for more 
details) . 

Using these assumptions it can be shown that, just as was the case for the structure function 
in Equation (j4]), the dispersion function 1 — (cos [A^ {£)]) can be decomposed into turbulent and 
ordered terms 



1 - (cos [A$ {£)]) = [62 (0) - 52 (£)] + (0) - {£)] 
= {b^{0) + [a^{0)-aH£)]}-b^i£), 



(13) 
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with the normalized ordered and turbulent autocorrelation functions given by 



{£) 



(Bq-Bq {£)) 
(B-B(0)> 

(BfBt(£)) 
(B-B(0)> 



(14) 
(15) 



respectively. The quantity 6^ (0) in Equation (jl3p is simply, from Equation (jlSp , the integrated 
turbulent to total magnetic energy ratio. It is also the equivalent of (A$^ (^)) (see Equation ([6])) 
when dealing with the angular structure function. The ordered function (0) — {£), which we 
assume to be of a larger spatial scale than 6^ (£), can be advantageously modeled with a Taylor 
series. Since, as we will soon discuss, we adopt a model of turbulence where the autocorrelation 
function is even in directions parallel and perpendicular to the projection of Bq on the plane of the 
sky, it follows that we can write 



(0) - {£) 



oo 

^ a2i,2j 
i+j=l 

i,j>0 



.p2ip2j 



(16) 



Accordingly, we will proceed by fitting the part within curly braces in Equation (I13p using 



6= (0) + [o^- (0) - (t)] = (0) + V a2,,2itf4' 



i+j = l 

i,j>0 



(17) 



to the data for high enough values of £ = \£\ where we expect 6^ {£) to be negligible (i.e., b"^ {£) 
dominates at lower values of £). We will then obtain 6^ {£) by subtracting the dispersion function 
data (i.e., the left hand side of Equation (jl3p ) from the aforementioned fit. This function is the 
equivalent to Equation ([7j) for x (^) defined in Section [2?T] for the angular structure function analysis. 

We now adopt a picture for rnagne tized turbulence consistent with current models, either in- 
comp ressible (jGoldreich &: Sridhaiill995l ) or compressible (ICho &: Lazarianll2003l : lKowal &: Lazarian 
That is, we will assume that some anisotropy is present in the magnetized turbulent au- 
tocorrelation function where it is expected that the correlation lengths in directions parallel and 
perpendicular to the ordered magnetic field are different. Although this is undoubtedly an idealiza- 
tion, we model the intrinsic autocorrelation function of the magnetized turbulence as an ellipsoid 
Gaussian function, with the symmetry axis of the ellipsoid aligned with the ordered magnetic field 
Bq (see Fig. [3]). It is possible to analytically solve for 6^ {£) for such cases with the further assump- 
tion that the telescope beam profile is circular Gaussian in form 
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Fig. 3. — Orientation of the Gaussian ellipsoid characterizing our idealization of magnetized turbu- 
lence in the {vi ,V2,u) coordinate system. The vi and V2 axes define the plane of the sky while the 
line of sight point along the negative ?i-axis. The inclination angle relative to the line of sight of 
the ellipsoid symmetry axis (and of Bq) is given by a while the angle /3 defines the orientation of 
its projection on the plane of the sky. 



with W = 0.425 FWHM the beam radius (i.e., it is not the beam's FWHM but its "standard 
deviation" equivalent). A general solution for 6^ {£) when the magnetized autocorrelation function 
ellipsoid is at arbitrary inclination relative to the line of sight and arbitrary projected orientation 
on the plane of the sky can be derived and is given in the Appendix. 

For M51 we will limit ourselves to three cases. One has the ellipsoid symmetry axis (and the 
ordered magnetic field Bq) at an inclination angle a relative to the line of sight and its projection on 
the plane of the sky advantageously aligned with one of the observers coordinate axes (still on the 
plane of the sky; see the Appendix). The integrated (or beam-broadened) autocorrelation function 
is then analytically expressed by 



-\[iy{5l+2w^)+ty{t4+2w^)] 



(19) 



N{Bl) + {B^,) 

where the number of turbulence cells contained in the gas probed by the telescope beam is given by 



{51 + 2W^) {fil + 2W^) A' 



(20) 



with 



5f, sin^ (a) + 5^ cos^ (a) 



(21) 



and S\\ and 6± the turbulent correlations lengths parallel and perpendicular to the (local) ordered 
magnetic field Bq, respectively. The inverse of these correlation lengths are, in effect, the cor- 
responding widths of the associated turbulent power spectrum (see the Equation (IA7p and the 
associated discussion in the Appendix). We once again emphasize that, in Equation (|19p . ii and £2 
are the displacements respectively perpendicular and parallel to the projection of Bq on the plane 
of the sky. In Equa tion (1201) A^ is the effective depth of the column of polarized gas probed by 
the telescope beam (jHoude et al.l l2009l : see their Section 3.2 and Equation (45)). This elliptical 
Gaussian turbulence case will actually be the last one considered in Section 14.2.21 



The first case considered (in Section 14. ip will be for isotropic turbulence when = 5j_ = 6; 
Equations (llOp and (|20l) then reduce to 
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and 



N = ^ =^ . 23 



Equation p6p is furthermore simplified to 



{Q) - {I) = Y,a2ji^^ . (24) 



This isotropic Gaussian turbulence is the model previously solved and used in lHoude et al.l (j2009l ) 



Also, in Section 14.2.11 we consider what could be termed as an "hybrid" model of anisotropic 
turbulence where Equations (122^ - (l24l) for isotropic turbulence are applied to independent analyses 
on displacement vectors that are grouped in sets oriented approximately parallel or perpendicular 
to the local mean magnetic field. Differences in the width of the two turbulent autocorrelation 
functions can then reveal anisotropy in the magnetized turbulence. 

For all three cases, the parameters measurable by fitting the integrated turbulent autocorrela- 
tion data with Equation (fTOl) (or Equation (122^ ) are the correlation lengths 5\\ and 5± (or simply 5), 
and the intrinsic turbulent-to-ordered magnetic energy ratio (-Bt ) / (-^o) ("^^ assume that A' and 
the inclination angle a are known a priori). 



3. Observations 



In this paper we use the high resolution radio polarization observations of lFletcher et al.l (j201ll ) . 
M51 was observed with the VLA at A6.2 cm using the C- and D-array configurations. Standard 
data reduction and imaging was carried out using AIPS to produce maps of the Stokes parame- 
ters I, Q, and U . These maps were merged with an Effelsberg 100-m telescope map at the same 
wavelength in order to correct for missing large-scale emission in the VLA data. The polarization 
angles we are using in this paper were calculated as <I> = ^ arctan (U/Q) and the polarized intensity 
as P = Y^Q^ + U'^, with a first-order correction for the positive bias in polarization due to noise. 
This is accomplished by simply subtracting the polarization uncertainty dp to the measured po- 
larization intensity Pm su ch that ~ P'^ — cTp, which is a good approximation when P > 3ap 
(jWardle fc Kronberg]|l974l ). We only use the highest resolution maps, with a FWHM of 4", a grid 



sampling of 1", and RMS noise of aj = 15 /uJ y/beam ^ in Stokes / and crp = 10 /iJy/beam ^ in P. 



At the assumed distance of 7.6 Mpc for M51 (jCiardullo et al.ll2002l ) 4" w 150 pc. The three regions 
shown in Figure [1] contain 520, 229, and 301 data points that verify P > 3ap, which are used for the 
corresponding analyses for the spiral arms in the northeast and southwest, and the galaxy center, 
respectively. 
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3.1. Data Analysis 



As stated in Sections 12.11 and I2.2| the polarization angle ^ is the basic o bservable needed fo r 
our analysis. Following the detailed discussion presented in Appendix B of iHoude et al.l (j2009l ). 
given the angle difference between a pair of data points separated by a distance 



ri - Ti 



(25) 



we calculate the (average) function (cos (A$jj))^ from the data for all (ik — A£/2) < iij < (l^ + A£/2), 
with £fc = kAi corresponding to an integer multiple of the grid spacing A£ = 1". This function is 
then corrected for measurement uncertainties according to 



(cos {A<^ij)) 



where the uncertainty on A^ij is given by 



(cos(A$ij-))fe 



fc,0 



l-i(a2(AcI>,,)),^ 



(26) 



aHA^^j) ^ ^\^^) + <7\<^,) - 2a(<I>,)^('J>,)e-^i/^^' 



(27) 



and (T^(<I>i) is the uncertainty on <l>j. Finally, the measurement uncertainties for the adopted dis- 
persion function 1 — (cos (A<&jj))^ ^ is determined with 



(cos (AcD,,)), J = (sin (A<I>,,))^ (^'(A^^,))^ + T (cos {A^^J))l (^t'(A<I>.,)>, , (28) 



for all (^fc — Ai/2) < lij < {1^ + Al/I). The data and results presented in the figures and tables 
that follow are all based on these equations. 



3.2. Faraday rotation 

Our analysis assumes that the observed polarization angles at A6 cm trace the orientation 
of the local magnetic field (we ignore the 7r/2 difference between the linear polarization plane of 
the observed electric field and the orientation of the magnetic field at the source of the emission). 
Faraday rotation can add an extra level of complexity to the distribution of angles, so here we 
estimate its contribution to the observed angles. 

Faraday rotation can produce a systematic variation of ^> with position due to the presence of 
a mean magnetic field; if the mean field lies in the same plane as the galaxy d isc then the pos i tiona l 
variation will occur due to the inclination of the galaxy to the line of sight. iFletcher et al.l (j201ll ) 
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modeled the mean magnetic field in M51 and found that it does lie in the galaxy plane and is weak. 
The maximum observed rotation measure due to the mean field of their model is RM ~ lOradm^, 
which rotates A6 cm emission by A$ ~ 2deg. 



Random fluctuations of Faraday rotation, ctrmj will also produce fluctuations in $. iFletcher et al 



(|201ll ) estimated that the intrinsic standard deviation of RM in M51 at 15" resolution is ctrm ~ 
lOradm^. At the 4" resolution we are using ctrm will be stronger, scaling as the ratio of the beam 
widths, so in our data ctrm ~ 40radm^ corresponding to a rotation of ~ 8deg at A6 cm. 

Thus Faraday rotation produces uncertainty in our dispersion functions of 1 — (cos(A<I>)) ~ 
0.01. This uncertainty is about an order of magnitude below the difference between the fits using 
Equation (I16p (or (I24p ) and the observations, until decorrelation of the angles occurs (e.g., see the 
middle panels of Figs. |l]to|6|. Note that decorrelation, i.e., where the autocorrelation function 
becomes zero, mostly occurs when the dispersion function is about 0.1, which corresponds to an 
angle difference of about 30deg. Therefore we will ignore the contribution of Faraday rotation to 
the polarization angles in our data, ot her than as a source of error. We also note that A3 cm data 



also presented in lFletcher et al.l (j201ll ). which suffer less from Faraday rotation, were not used for 



this analysis because of lower spatial resolution. 



4. Results 

In this section we present the results of the angular dispersion analyses for isotropic and 
anisotropic magnetized turbulence. We used the polarization data from the three regions iden- 
tified in Figure [U the northeast spiral arm, the centre of the galaxy, and the southwestern spiral 
arm. In all cases, we only consider measurements for which p > Scjp, where p = P/I and Gp are the 
polarization level and its uncertainty, respectively. 



4.1. Isotropic Turbulence 

We first consider the case of isotropic turbulence and model our data for the three suitable 
regions in M51 with Equations (|22l) - (i24| l. All the pertinent functions are therefore assumed to 
possess an azimuthal symmetry about the £ = axis. 

Figure U] shows the result of the isotropic dispersion analysis for the northeast spiral arm as a 
function of £^ (top) and I (middle). The broken curve ("ordered") is the least-squares fit for the 
sum of the integrated turbulent-to-total magnetic energy ratio (6^ (0) in Equation ([22])) and the 
ordered component of Equation (j24p (i.e.. Equation (|17p while using Equation (|24p on the right- 
hand side) to the data contained within 6 < £ < 10; data points are shown with symbols. The 
integrated magnetized turbulence autocorrelation function iP' {(), obtained by subtracting the data 
from the aforementioned fit of the middle graph, is shown at the bottom. The broken curve on the 



- 16 - 



bottom graph shows the radial profile of the "autocorrelated synthesized beam." This represents 
the contribution of the synthesized beam to the (width of) 6^ (£). That is, this is what 1? {£) would 
look like in the limit where the intrinsic turbulent correlation length 6 were zero in the exponent of 
Equation (p2|) (i.e., disregarding its effect on the amplitude of (i) through N). It follows from this 
and the fact that the data points for 6^ {£) fall practically on top of the autocorrelated beam that the 
correlation length 6 in this region of M51 is significantly smaller than the beam size W ~ 1'.'70 ~ 63 
pc. We also find from these graphs that 6^ (0) ~ 0.028, however, we cannot proceed any further in 
view of the impossibility of determining 5 for this data set. 

Figures [5] and [6] show the results of the dispersion analyses for the centre and the southwest spiral 
arm of M51, respectively. In both cases we can clearly see a broadening of the integrated magnetized 
turbulence autocorrelation function beyond that due to the telescope beam (bottom graphs); this is 
an imprint of the turbulent correlation length intrinsic to the magnetized turbulence. Least-squares- 
fitting a Gaussian function to these reveals that 5 = 2'.'0 (74 pc) and I'.'T (62 pc), respectively. It 
therefore follows that we can provide estimates for the number of turbulent cells probed by the 
telescope beam and the intrinsic turbulent-to-ordered magnetic energy ratio for these t wo re gions. 
The results are presented in Table [H where we set A' = 80 pc from iFletcher et al.l (j201ll ) . We 
thus find that our results are in good agreement with those of lFletcher et al.l (j201ll ) who estimated 
Bt/Bf) ~ 1 in the neighborhood of the spiral arms using a rotation measure dis persion analysis. 
Similar values have been reported in previous analyses for other sources (e.g., see iBeck et al.lll9l 
for NGC 6964). As will be discussed in Section [5l their value of that 25 ~ 50 pc is consistent to 
ours given the uncertainty in some of the parameters that enter the analysis. 



4.2. Anisotropic Turbulence 

4-2.1. "Hybrid" Model of Anisotropic Turbulence 

We now abandon the isotropy assumption and we make a first attempt at treating the more 
general case of anisotropic turbulence. To do so, we define two separate bins of data where the 
polarization angle differences used in Equation (jlOp are such that the displacement i is either 
oriented within ±45° of the (plane of the sky component of the) local mean magnetic field (and 
labeled .^n) or in a direction within ±45° from the axis normal to it (labeled £±^). This is illustrated in 
Figure [71 The orientation of the mean field at a given position is simply approximated by averaging 
polarization angles contained within a radius of 2". We then perform two separate dispersion 
analyses that will allow us to measure differences in the magnetized turbulence correlation lengths 
parallel and perpendicular to the local field, and 5±, respectively. Although this analysis is not 
based on the more rigorous model given in Equations (llOp and (i20l) . it will allow us to look for direct 
evidence of anisotropic turbulence in the same three regions of M51, as was done in the previous 
subsection. A more rigorous analysis applied to M51 as a whole will follow in Section [4.2.21 

Figure [8] shows the results for the northeast spiral arm previously analyzed under the isotropy 
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Fig. 4. — Isotropic dispersion function for the northeast spiral arm as a function of (top) I and 
(middle). The broken curve ("ordered") is the least-squares fit for the sum of the turbulent to 
total magnetic energy ratio and the ordered component to data contained within 6 < £ < 10; data 
are represented with symbols. Bottom: the magnetized turbulence autocorrelation function 6^ (£) 
obtained by subtracting the data from the aforementioned fit of the middle graph. The broken 
curve shows the radial profile of the "mean autocorrelated synthesized beam." 
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Fig. 5. — Same as Figure |4] but for the centre of M51. 
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Fig. 6. — Same as Figure |4] but for the southwest spiral arm of M51. 
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Table 1. Results for isotropic turbulence. 





Northeast Arm 


Centre 


Southwest Arm 


6 (pc)'^ 




67 ± 7 


66 ±8 






13 ±3 


14 ±4 


{Bi)/{By 


0.028 ± 0.002 


0.088 ± 0.026 


0.072 ± 0.025 


{Bl)/{BlY 




1.28 ±0.29 


1.08 ±0.29 






1.13 ±0.13 


1.04 ±0.14 



^Turbulent correlation length (1" = 37 pc); from the fit of Equa- 
tion (j22p to the data. 

^Number of turbulent cells probed by the telescope beam, using 
A' = 800 pc; from Equation (f23l) . 

*^Measured value for the integrated turbulent to total magnetic 
energy ratio, corresponding to b"^ {£ = 0) = (B^) / [N (-Bq) ± (Bt)] 
(see Equation (|22]l ). 

'^Turbulent to ordered magnetic energy ratio, corrected for signal 
integration; from the fit of Equation ([22|) to the data. 

'^Calculated from the root of (S^) / (B^). 
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Fig. 7. — Grouping of differences A<I> {£) into two sets depending whether or not £ is oriented in 
a cone whose boundaries are within ±45° of the orientation of the local mean magnetic field. The 
displacement vectors are labeled with i± or £|| depending on the case (i.e., perpendicular or parallel 
to the field, respectively). 



- 21 - 



assumption in Figure [H For such analysis the dispersions function parahel and perpendicular to 
the magnetic field must be treated simultaneously. That is, the least-squares fits for the sum of 
the integrated turbulent-to-total magnetic energy ratio and the ordered component (broken curves 
in the top two graphs) to the data contained within 6 < ^ < 10 are not independent since they 
must meet at £|| = = 0. These fits are thus performed simultaneously. The bottom graph 
shows the integrated magnetized turbulence autocorrelation functions parallel and perpendicular 
to the mean field as well as the mean autocorrelated telescope beam, as before. Although we see 
evidence for anisotropy from the separation of the two autocorrelation functions, we find that the 
perpendicular function has a width that is narrower than the contribution of the telescope beam, 
which is impossible. This is most likely due to the fact that our fits (on the top two graphs) are made 
with data points that are located at too low values for the displacements £|| and l± and therefore to 
some extent fail to cleanly separate the ordered and turbulent dispersion functions (at the expense 
of the latter; see Sec. [5|). At any rate, we can infer from this analysis that ^-^t ^ / {^0^ ~ ^'^"^ (™ 
agreement with the isotropic analysis) and that the intrinsic magnetized turbulence autocorrelation 
appears to be broader along the local magnetic field orientation than perpendicular to it. 

Figures [9] and [10] show the results of the same anisotropic dispersion analysis for the center 
and southwest spiral arm of M51, respectively. In these two cases, however, we clearly resolve the 
anisotropy in the turbulence. That is, we observe a separation in the integrated autocorrelations 
functions (bottom graphs) along the directions parallel and perpendicular to the local mean magnetic 
field. The former being the broader of the two, which is also consistent with was observed for the 
northeast spiral arm in Figure [HI As we will discuss in Section [5| this result is consistent with th eory 



and simulations of incompressib le (iGoldreich Sridhai 



1995; 



Cho. Lazarian. &: Vishniadl2002l) and 



compressible MHD turbulence ()Cho &: LazarianI [20031 : iKowal fc Lazarian 



201 



The level of anisotropy in the turbulence can be gauged through the parallel-to-perpendicular 
correlation length ratio 5^^/6±, which is measured to be approximately 1.3 and 1.5 for the centre and 
southwest spiral arm, respectively. To estimate the number of turbulent cells N we use Equations 
(|20p and (|2ip . with a = tt/2. We once again find that the turbulent-to-ordered magnetic field 
strength ratio is significant and hovers around unity with 0.8 < Bi/Bq < 1.7. A summary of the 
results is given in Table [2j 



4-2.2. Two-dimensional, Anisotropic Gaussian Turbulence 

We finally perform one last anisotropic analysis by taking advantage of the large number of 
reliable polarization measurements contained in the complete map of M51 shown in Figure [TJ That 
is, we now consider the whole map at once without discriminating between the different regions 
(as long as p > 3ap). We hope in doing so that the large number of measurements will allow for 
the characterization of the intrinsic two-dimensional turbulence autocorrelation function, using the 
Gaussian anisotropic model given in Equations (llOp and (I20p . One would expect that the previously 
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Fig. 8. — Same as Figure|4]for the northeast spiral arm, but for directions parahel and perpendicular 
to the local mean magnetic field. 
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Fig. 9. — Same as Figure [8] but for the center of M51. 
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Fig. 10. — Same as Figure [8] but for the southwest spiral arm. 
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Table 2. Results for anisotropic turbulence. 





Northeast Arm 


Centre 


Southwest Arm 




25 ±9 


111 ±7 


61 ± 7 






87 ±8 


41 ± 5 






1.27 ±0.14 


1.48 ±0.25 


Arc 




7± 1 


32 ±8 


{Bi)/{By 


0.028 ± 0.002 


0.093 ± 0.003 


0.082 ± 0.001 


{Bl)/{BlY 




0.68 ±0.10 


2.86 ±0.68 


Bt/Bo' 




0.83 ±0.06 


1.69 ±0.20 



^Turbulent correlation length parallel to Bq (1" = 37 pc); from 
the fit of Equation (f22|) to the data. 

^Same as note a, but perpendicular to Bg. 

'^Number of turbulent cells probed by the telescope beam, using 
A' = 800 pc; from Equation (l20l) . 

'^Measured value for the integrated turbulent to total magnetic 
energy ratio, corresponding to b'^ {i = 0) = (B^) / [N (Sq) ± (Bt)] 
(see Equation (jl9p ). 

^Turbulent to ordered magnetic energy ratio, corrected for signal 
integration; from the fit of Equation (j22p to the data with S = 

^Calculated from the root of (B^) / (B^). 
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measured anisotropy, quantified with S^^/6±, would become more pronounced since we would do away 
with the cone-averages exemplified in Figure [71 

Figure [11] shows the result of the two-dimensional dispersion analysis. Only one quadrant 
of the contour plot of the dispersion function is displayed since it is assumed even in directions 
perpendicular and parallel to the mean magnetic field (i.e., in even powers of ti and £.2)- Since we 
must now least-squares fit a two-dimensional surface corresponding to the sum of the turbulent-to- 
ordered magnetic energy ratio and the ordered component of the dispersion function (i.e., the right- 
hand side of Equation (jl3p ). it is to be expected that this fitting process will be more challenging 
than before. This can be verified in Figure [12] were cuts through the two-dimensional dispersion 
and integrated turbulence autocorrelation functions along directions parallel and perpendicular 
to the local mean magnetic field are shown; data contained within 7 < ^ < 10 were used to 
perform the aforementioned least-squares fit to the right-hand side of Equation ()13p . A comparison 
of Figure [12] with any such figures stemming from the previous isotropic or anisotropic analyses 
reveals that our fit to the two-dimensional dispersion function (top and middle graphs) is unable to 
perfectly match the data. The main consequence of this being the somewhat "ragged" appearance of 
the integrated two-dimensional turbulence autocorrelation function presented in the bottom graph 
of Figure [12] (symbols). Nonetheless, it is interesting to note that we once again find the same 
anisotropy as before, with the result that (Jy > 5±. We sought to quantify this by performing a two- 
dimensional (elliptical) Gaussian least-squares fit, using Equation (|19p . to the integrated turbulence 
autocorrelation function data. This is shown in the top panel of Figure [T3l where a contour plot 
of the aforementioned Gaussian fit (red) is superposed on that of the integrated two-dimensional 
turbulence a utocorrelation function (black). For this we used the known value of a = 70° for M51 
(lTullylll974l ). The fit is forced to be even in directions perpendicular and parallel to the local mean 
magnetic field (i.e., even in powers of ii and £2; see Equation (flBl) ). as the dispersion function 
was also assumed to be. Although this Gaussian fit appears to be reasonable for i < 4", it is 
not expected to be a realistic representation since it is unlikely that the magnetized turbulence 
is Gaussian in nature in M51. Nonetheless, it allows us to extract a useful approximation to the 
intrinsic two-dimensional magnetized turbulence autocorrelation function; the resulting function is 



shown in the bottom panel of Fig ure [T3l As stated ear . 



and simulations of incompressibl e (iGoldreich &: Sridhar 



compressible MHD turbulence (ICho &: Lazarian 



2003; 



ier, such results are consistent with th eory 



1995; 



Cho. Lazarian. &: Vishniadl2002n and 



Kowal fc LazarianI I2OIOI ) . The parameters 



extracted from this anisotropic analysis are also consistent with our previous results for the three 
independent regions and are presented in Table [3] 

Our results do not take into account any systematic uncertainties on some of the parameters 
used to characterize M51. For example, the effective depth A' = 800 pc, which comes in for all 
three cases treated in this section (isotropic, "hybrid" anisotropic, and anisotropic turbulence), enters 
linearly in the evaluation of the number of turbulent cells contained in a telescope beam (see 
Equations (I20p and ((23l) ). In turn, the relative strength of the turbulent magnetic field component 
to the ordered magnetic field scales inversely with N^/"^ . An overestimation by a factor of two in 
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A', for example, would bring a corresponding underestimate of Bt/Bo by \/2. Furthermor e, the 
fully a nisotropic model is also dependent on the inclination of the galaxy, which according to iTullv 
(jl974l ) spans a = 70° it 5°. Unlike its dependency on A' discussed above, the relative level of 
turbulence is found to be largely insensitive to changes in a. On the other hand the correlation 
lengths are somewhat affected by such uncertainties. For example, we find that 100 pc > > 96 pc 
and 1.92 > £\\/e± > 1.78 when 65° < a < 75°. 



5. Discussion 



Our application of the dispersion analysis of iHildebrand et al. ( 2009 ) and Houde et al. ( 2009 



2011 



, and its generalization to include anisotropy, to M51 reveals some interesting information on 
the nature of magnetized turbulence in this galaxy. As was prev iously mentioned, it is important 
to note that both our analysis and that of iFletcher et al.l (j201ll ) yield results that are consistent 



with one another, ev en though they are based on completely different approaches. For example. 



Fletcher et al.l (120111 ) determined the size of a turbulent cell (i.e., approximately twice the turbulent 



correlation length) by using their measured dispersion of rotation measure, while accounting for 
the averaging of turbulence inherent to the observation process (see their Equations (3) and (5)). 
Their value of approximately 50 pc for the size of a turbulent cell can be readily compared with 
our estimates of (5 ~ 65 pc determined for isotropic turbulence in Section 14.11 (see Table [T]) . It is 
noteworthy that both techniques provide results that are within a factor of two or so from each 
other, which is interesting considering the uncertainty in the adopted values for some parameters 
(e.g., the mean electron density {rif.) =0.1 cm~^ used in their calculations). 



Fletcher et al.l (120111 ) were also able to discern between the contribution of the different com- 
ponents of the total magnetic field. They found that the total magnetic field i? ~ 15 /iG is split into 
an ordered (Bq) and an isotropic (i.e., random, Bt) components each of ~ 10 /xG. It is important 
to note that their definition of an "ordered" magnetic field is not restricted to a "mean" field, which 
would result from the average of the magnetic field vector over some suitable (large) scale. More 
precisely, they define the ordered field as that which is traced by the polarized emission. For M51 
they find that the ordered magnetic field is composed of a weak ~ 2 /iG mean component and an 
anisotropic random field of ~ 10 fiG (we note that such a "mean component" implies a field with 
a coherent direction, while the anisotropic field has reversing directions). This anisotropic field, 
presumably resulting from "compression in the spiral arms or localized enhanced shear," would then 
display a stro nger azimuthal variation and thus be responsible for most of the polarized radio emis- 



sion in M51 (IFletcher et al.ll201ll ). Their analysis therefore yields Bi/Bq ~ 1, as was mentioned 
earlier. Since our dispersion analysis is based on changes in the orientation of polarization vectors 
with length scale, it will not be able to discriminate between anisotropic random and mean field 
components of the ordered field and we cannot comment on its detailed nature. We note, however, 
that our estimates B^/Bq ~ 1.13 and 1.04 determined fo r the i sotropic turbulence case in Section 



14.11 are in excellent agreement with that of lFletcher et al.l (|201ll ) (see Table [T]). The typical degrees 
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Fig. 11. — Contour plot of the two-dimensional dispersion function for the whole M51 polarization 
map of [TJ This function is assumed to be even in directions parallel or perpendicular to the mean 
magnetic field. 
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Fig. 12. — Top and middle: cuts through the two-dimensional dispersion function (symbols) and the 
ordered fits (broken curves; using values of 7 < ^ < 10) along directions parallel and perpendicular 
to the mean magnetic field. Bottom: the corresponding profiles for the turbulence autocorrelation 
function. 
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Fig. 13. — Top: contour plot of the two-dimensional turbulence autocorrelation function (black) 
and a Gaussian fit (red); the fit is forced to be even in directions parallel and perpendicular to 
the mean magnetic field. The contours are drawn at 10 to 90 percent (10 percent increments) of 
the peak iP' (0) = 0.063. Bottom: the intrinsic turbulence autocorrelation function, where we set 
W = a = 7r/2, and = 1 in Equations (|19p and (|2ip . while using the values for 5\\ and 5± 
obtained with the Gaussian fit (red) in the plot on the top graph (see Table [3]). The contours are 
drawn at 10 to 90 percent (10 percent increments) of the peak (-B^) / (-B^) = 0.50. 
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Table 3. Results for two-dimensional anisotropic turbulence (a = 70°). 



All regions 


5|| (pc)^ 


98 ±5 


S± (pc)^ 


54 ±3 




1.83 ±0.13 




15 ±2 


{Bi)/{By 


0.063 ± 0.008 




1.02 ±0.08 


Bt/Bo' 


1.01 ±0.04 



^Turbulent correlation length 
parallel to Bq (1" = 37 pc); 
from the fit of Equation (|19p to 
the data with a = 70°. 

^Same as note a, but perpen- 
dicular to Bq. 

'^Number of turbulent cells 
probed by the telescope beam, 
using A' = 800 pc; from Equa- 
tion (|20ll . 

'^Measured integrated 
turbulent to total mag- 
netic energy ratio, corre- 
sponding to 6^ = 0) = 
(B?) / [N {Bi) + {B^)] (see 
Equation (fTOlll. 

''Turbulent to ordered mag- 
netic energy ratio, corrected for 
signal integration; from the fit 
of Equation (llOp to the data 
with a = 70°. 

^Calculated from the root of 
{B?)/{Bi). 
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of polarization, at A6 cm and 4" resolution, in the three regions that we consider are p = 34% in 
the northeast arm, p = 29% in the southwest arm, and p = 15% in the centre: in ca lculating these 



values we have assumed that 25% of the continuum emission at A6 cm is thermal (Fletcher et al. 



20 111 ). For the arms these data are in good agreement with the degree of polarization expected when 
Bt/Bo « 1, as p = 0.7 B^/{B^ + Bf) 0.35. This shows that the estimates of Bt/Bo derived from 
our polarization angle dispersion analysis are compatible with independent methods of interpreting 
polarization data. In the centre of M51 the physical environment is different from the rest of the 
disc, due to the presence of an AGN and jet, and so the fraction of thermal emission and thus the 
degree of polarization are harder to estimate. In addition, in this region the polarized emission can 
originate from a different location than a large fraction of the total synchrotron emission, again 
making a useful estimate of p difficult. 

Perhaps the most important result stemming from our analysis is the clear detection of ani- 
sotropy in the magnetized turbulence. Whether we consider the three analyzable regions of M51 
separately or together we consistently find that the turbulent correlation length is larger in a di- 
rection parallel to the mean orientation of the local magnetic field than in a direction perpen- 
dicular to it (see Tables [2] and [3]). As was mentioned in Section 14.21 this result is predicted by 



2002 



Cho &: Lazarian 



2003 



Kowal &: Lazarianll201Cl ). Such anisotropy 



current theories for incompressible and compressible MHD turbulence (IGoldreich &: Sridhai 
Cho. Lazarian. fc Vishniad 



1995 



has a lso been observed in the Taurus (jHever et al.ll2008l ) and Orion molecular clouds (jChitsazzadeh et al 
2OI2I ) within our Galaxy. The level of anisotropy we observe in M51, which we quantify with the 
parallel-to-perpendicular correlation length ratio d\\/S±, is significant and in qualitative agreement 
with numerical simulations of magnetized turbulence. Our contour plot of the turbulent autocor- 
relation function (see the bottom panel of Figure [T3l) can be compared with the simulations of 
incompressible MHD turbulence of ICho. Lazarian. &: Vishniad ([2002), for example. In particular, 
the turbulent velocity correlation function presented in their Figure 6 has an appearance that is 
similar to our derived intrinsic two-dimensional magnetized turbulence autocorrelation function. 
We would expect such similarities between these two types of autocorrelation functions under the 
flux-freezing approximation, which should hold in the medium probed with synchrotron polariza- 
tion observations. Our clearest measure of anisotropy uses the more comprehensive Gaussian model 
defined with Equations (jl9p to ([2T]) . where the polarization data were analyzed for M51 as a whole 
without discriminating between the different regions. The significant amount of anisotropy thus 
measured, with — 1.8, is a statement of the importance of magnetic fields on the dynamics 

of the gas probed by the observations. We should also note that this anisotropy of t h e mag netized 
turbulent autocorrelation function is different from that discussed by lFletcher et al.l (120111 ). which 
pertains to the relative intensity of the two orthogonal components of the random magnetic field. 
Our analysis cannot say anything concerning any such anisotropy in the field strength, it can only 
inform us on the relative turbulent energy contained in the magnetic field through measurements 
of (B^) I (Bq), for example. 
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5.1. Shortcomings of the Dispersion Analysis 

Although the quahty of the data and the high resolution with which they were obtained allowed 
us to determine some fundamental parameters that characterize magnetized turbulence in M51, we 
expect that a slightly higher resolution and samp ling rate would result in an even more exhaustive 



analysis. As was shown by iHoude et al.l (l201ll ) using submillimeter dust polarization data for 
Galactic molecular clouds, spatial resolutions resulting in smaller telescope beams such that 5 > 
\plW not only allow the determination of the same parameters uncovered by the present analysis, 
they can also potentially reveal the underlying turbulent power spectrum. This is because the 
beam-broadened turbulent autocorrelation function {€) is related to the turbulent power spectrum 
\? (kt,) through a simple Fourier transform. It is then found that 



B 



^ y 7^t (k„, ku) sinc^ {~^^ 



(29) 



with TZi (k„, ku) and H (k^) the Fourier transforms of the intrinsic turbulent autocorrela tion function 



(i.e., not beam-broadened) and the telescope beam, respectively (see the Appendix, and lHoude et al. 



2011 



for a detailed discussion). It follows that beams of smaller spatial extent than the intrinsic 
turbulent autocorrelation function will have a broader spectral coverage that will reduce their fil- 
tering effect on the power spectrum. It then becomes possible to effectively invert Equation (|29|) to 
reveal the underlying power spectrum (through some "deconvolution" techniques, for example). In 
such cases, it is not necessary to assume any model for the turbulence, such as the Gaussian form 
used in our analysis. The measured turbulent power spectrum could thus be modeled directly from 
the data and compared to candidate theories for magnetized turbulence. 

As can be seen from our results for 6^ {£) in Figures Bll6l [8lfT0| and ll2l however, the contribution 
of the correlation length to the width of the beam-broadened turbulent autocorrelation function 
(approximately gauged through the ratio 5^/ (5^ -|- 2VF^); see Equations (|20p and (|23p ) is typically 
modest, implying that the spectral filtering of the telescope beam is too severe to recover the intrinsic 
turbulent power spectrum. But even a relatively modest increase in spatial resolution, e.g., by a 
factor of a few, could allows us to recover the power spectrum in future observations. 

Another negative impact of a larger telescope beam and its broadening of the autocorrelation 
function {£) is that it renders more difficult the separation of the small and large scale components 
present in the dispersion function. For M51 this means that the scale of the turbulence, quantified 
with 6, can get "mixed up" with the larger scale of the spiral structure through its artificial broad- 
ening to 6'^ + 2W'^ caused by the beam. As alluded to in Section 14. H this may be a reason why we 
were unable to see any contribution from S to the width of 6^ {£) in our analyses of the northeast 
arm (see Figs. [4] and [8]). More precisely, we were unable to cleanly separate the large from the small 
scale using our Taylor expansions, i.e., Equations (iT6|) and (I24p . This is probably also true, but to a 
lesser extent, for the other two regions studied, as can be seen from the absence of significant "tails" 
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for I > 5" in 6^ (£). An increase in spatial resolution would resolve this issue, which is likely to 
bring some error in our determination of the correlation length scales and turbulent to total energy 
ratios. We expect this error to be small, but it is not possible to quantify it at this point. 

Finally, we wish to once again emphasize that the choice of a Gaussian turbulence model is 
unlikely to be realistic for M51. But in view of the aforementioned impossibility to uncover the 
underlying turbulent power spectrum because of the significant spectral beam filtering, this model, 
which can be solved analytically, allows us to quantify key parameters that characterize magnetized 
turbulence. Furthermore, our more comprehensive model for anisotropic Gaussian turbulence de- 
fined with Equations (I19p and (j20p implicitly assumes that the N ellipsoid turbulent cells contained 
in the column of gas probed by the telescope beam have the same spatial orientation in relation to 
the local magnetic field. This is clearly unlikely to be true across the beam (FWHM ~ 148 pc), or 
through the thickness (~ 800 pc) and the extent of the studied regions on the galactic disk. It is 
therefore more realistic to view the correlation lengths and 5± as some averages representative 
for magnetized turbulence in M51. 



6. Summary 



We conducted a dispersion analysis, using a generalization of the tech nique of 



([20091) to previously published high-resolution synchrotron polarization data (jFletcher et al 



Houde et al. 



2011 



with the goal of charactering magnetized turbulence in M51. We first analyzed three distinct regions 
(the center of the galaxy, and the northwest and southwest spiral arms) and measured the turbulent 
correlation length scale from the width of the magnetized turbulent correlation function for two 
regions and detected the imprint of anisotropy in the turbulence for all three. Furthermore, analyzing 
the galaxy as a whole allowed us to determine a two-dimensional Gaussian model for the magnetized 
turbulence in M51. We measured the turbulent correlation scales along and perpendicular to the 
local mean magnetic field to be, respectively, 5\\ = 98ib5 pc and 5±_ = 54zb3 pc, while the turbulent 
to ordered magnetic fie ld strength rat i o is fo und to be B^/Bq = 1.01 it 0.04. These results are in 
agreement with those of lFletcher et al.l (|201ll ). who performed a Faraday rotation dispersion analysis 
of the same data. Finally, our detection of anisotropy, quantified with a parallel-to-perpendicular 
correlation length ratio with 5\\/5±_ ~ 1.83 it 0.13, is consistent with current magnetized turbulence 
theories. 
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A. Anisotropic Gaussian Magnetized Turbulence 

The cloud- and beam-integrated magnetic field is defined with 
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B(r) 



jj H{T-a) i ^ Fo (a, B (a, z) dz 



(Al) 



where the beam profile is denoted by H (r), while the weighting function i<o (a, z) > is the 
(ordered) polarized emission associated with the magnetic field B(a, z), and A is the maximum 
depth of the cloud along any line of sight. The integrated autocorrelation function is then 



(BfBt(£)) 



H (a) H (a' + i) 



u 



1 - — ) 7?.t (v, u) du 



d^a'd^a, 



(A2) 



with nt(v,u) = (Fn(a, z) Fq (a', z')) (Bt (a, z) • Bt (a', z')), u 



and V 



1)262; iHoude et al.ll2009l ). We refer to {vi,V2-,u) as the observer coordinate system, the vi and V2 
axes define the plane of the sky, while the line of sight point along the negative u-axis. 

The assumptions of statistical independence between the turbulent and ordered magnetic fields, 
homogeneity in their strength across the source, and of overall stationarity previously stated in 
Section 12.21 are all required to arrive at Equations ()Aip and ()A2p . Of these, the assumption of 
homogeneity is particularly useful for analyzing our data. This is because synchrotron polarization 
signals bring in the complication that the weighting function Fq (a, z) is also a function of the 
magnetic field strength (approximately proportional to its second power), and would therefore 
appear to significantly jeopardize any calculations stemming form Equation ()A2p . However, this 
dependency is seen to disappear in the calculation of the angular dispersion function (Equation (jlOp ) 
when homogeneity is assumed, since this weighting function will have the same proportionality factor 
(due to the field strength) at all points in the source (i.e., in the integrands of Equations (lAip and 
()A2p ). Our analysis can then proceed in a manner similar to the simpler case of polarization dust 
emission signal s, where there is n o linked between the value of Fq (a, z) and the strength of the 
magnetic field (jHoude et al.l 120091 ). 



In cases where the autocorrelation function for anisotropic magnetized turbulence is idealized 
with a Gaussian ellipsoid we write 




where ^ (= (^1,1^2)) is a two-dimensional displacement vector perpendicular to the orientation 
of the local ordered magnetic field Bq and r is the displacement along Bq; a prolate example 
is shown in Figure [3] along with the relationship between the {vi,V2-,u) and {ii-,i2-,T) coordinate 
systems. We adopt a model for anisotropic magnetized turbulence where the symmetry axis of 
the ellipsoid is aligned with Bq; this also implies that the length scale of the ordered field is much 



- 36 - 



larger than the correlation lengths and 6± characterizing the turbulent field Bt. The function 
i^Q (^,r)) is the autocorrelation of the ordered polarized emission, which we approximate to a 
constant (-Fq) ^ (-^0 (O!*-*)) Equation ()A3P as it is assumed that its correlation length is also 
significantly larger than and 5±. We seek to express this function (i.e., Equation (jASP ) using 
the observer coordinates {vi,V2,u). Referring to Figure [3l the inclination angle relative to the line 
of sight of the ellipsoid symmetry axis (and of Bq) is given by a, while the angle /3 defines the 
orientation of its projection on the plane of the sky. The precise relationship between the two 
coordinate systems is 



such that 



with 



^1 = ui cos (/3) + U2 sin (/3) 

^2 = —vi cos (a) sin (/3) + V2 cos (a) cos (/?) + u sin (a) 
T = vi sin (a) sin (/3) — V2 sin (a) cos (/3) + u cos (a) , 



(A4) 



I U ui U2 ^ ^ ^ 
-:o+Tt- = ^ + ^ + ^ ^'VYV2 + -^UV\ ^ 



"12 



tUV2, 



(A5) 



cos^ (a) sin^ (a) 

To ~l~ 



err 



^ - ^ j sin (a) cos (a) sin (/?) 



(7? 



— \ sin (a) cos (a) cos (/3) 

On ' 



sin^ (a) sin^ (/?) cos^ (/3) + cos^ (a) sin^ (/3) 



+ 



(A6) 



sin^ (a) cos^ (/3) sin^ (/3) + cos^ (a) cos^ (/3) 



^2 



+ 



"12 



^ - ^ I sin^ (a) sin (/3) cos (/3) . 



Inserting Equations (jASp and (|A6p in Equation (jA3p we can express the turbulent autocorrelation 
function with a dependency on (vi ,V2,u). 
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It is advantageo us to solve Equation (IA2 ) by considering its Fourier transform (i.e., the turbu- 



lent power spectrum; iHoude et al. 



20091 . l201lh 



BtiK) = \\H{K) 



1 

2^ 



Tit (k„, ku) sinc^ ' 



dkii 



(A7) 



with the correspondence (v,u) ^ (kjj,ku) between the two domains, and then recover the autocor- 
relation function through the inverse Fourier transform 



(27r)2 // 



(A8) 



We note that the two-dimensional power spectrum given by Equation (IA7[) could be compared to 
a Kolmogorov-type spectrum, for example, by fir st multiplying it by 27rfc^ , (for a three-dimensional 
spectrum a factor of 47r/c^ would be required; see lHoude et al.ll2009l . l201ll ). 



The solution for this problem mostly rests on the repeated application of the following relation 
for the integration of Gaussian functions 



a 



{b/2aY 



(A9) 



with a and b some constants. Considering Equation ()A9|) and the fact that the Fourier transform 
of a rotated function equals the rotated version of the Fourier transf orm of the unrotated function 
(i.e., when a = /3 = 0; see Appendix 4 of lHoude &: VaillancourtI 120071 ) . the Fourier transform of the 
turbulent autocorrelation function can be calculated to be 



T^t (kt,, ky) 



-\{e^kl+^J.\kl+^llkl-2^ll2klk2+2-|lkukl-2-^lkuk2) 



(AlO) 



with 



e2 


= 5| cos 


ll 


^ (4- 


ll 


^ (4- 


A 


= (5j sin^ 




= (5jsin^ 


1A2 





6?, cos^ (q) -|- 5'^ sin^ (a) 



sin (a) cos (a) sin (/3) 



51 



sin^ (a) sin (/3) cos (/3) 



(All) 
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and k„ = kiei + ^262- 

We also note that because the depth of integration along the line of sight is expected to be 
much larger then the turbulent correlation lengths (i.e., A ^ 5\\ and A ^ 6±) we have 



e 2 



Sine 



2 ( '^u 



sine 



2 ( 



A:,, A 



27r 
A^' 



(A12) 



Inserting Equation ()A12p in Equation ()A7p with 



(A13) 



and further using Equation ()A9p . it is found that 



A 



Calculating the inverse Fourier transform of Equation (IA14j) . still using Equation (IA9j) . then yields 



(BfBt(€)) 



AT/ 



(A15) 



with 



9 {£;^\\,S±;a,P) 



Ull (^i + 2W^2) (A16) 



and 



1 



27r5||52 



{\A + /^l2/ (^i + 2^2) + 2Ty2] (^2 + 2H^2) A 



(A17) 



Following the treatment of Houde et al. (j2009l ) for the ordered component of the autocorrelati 
function we write 



ion 



(Al^ 
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where A' < A is the effective depth over which the signal is integrated along the line of sight (in 
our case approximately the thickness of the disk of M51), wh ich is closely relate d to the correlation 
length of the ordered polarized flux (see Sees. 2.3 and 3.2 in lHoude et al.ll2009l ). 

Combining Equations (lAlSP - IIAlSP we find 



(BfBt {£)) 
(B-B(0)> 



(A19) 



with the number of turbulent cells N probed by the telescope beam given by 



1 

N 



(A20) 



If we choose to align one of the axes of the observer's coordinate system on the plane of the 
sky (i.e., the W2-axis) with the large scale magnetic field Bq, then /3 = and Equations (jAlOP and 
(|A20p reduce to Equations (|19p and ()20p . respectively. 
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